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Abstract 



■ We present an explicit OPE analysis for the first moment of gi up to order 

(y-j I M^/Q^. This result allows to calculate power corrections to the Bjorken and EUis- 

00 ■ Jaffe sum rules. 

. 

PACS numbers: 13.88. +e, 12.38.Bx, 12.38.Qk 

m 

, The investigation of the spin structure of the nucleon is one of the most active helds 

>-H ' in hadron physics. This is demonstrated by the fact that CERN, SLAG, and DESY ah 

have programs in polarized lepton-hadron scattering. So far data are available from 
gn! EMC (proton) Q, SMC (deuteron) §, and E142 (helium) g. As ah of these groups 

(— I I took data at different values of these results can not be compared directly. In partic- 

ular for E142 the lowest x points are measured at only ^ 1 GeV'^ such that higher 
twist corrections should be important. Actually one has to face a typical track-off. At 
low one can get very good statistics but the theoretical analysis is more involved, 
' at high theory is easy, but the statistics of the experiment does not allow detailed 

statements. 

Polarized scattering has, however, the great advantage that the higher twist corrections 
seem to be calculable with a sufficient precission As also the HERMES experiment 
at HERA will run at small it is a necessary task for theory to clarify the issue of 
dependence. The most basic block for this is the operator product expansion for the 
moments of gi and §2- Unfortunately there was a mistake in the original work Q as 
noticed by Ji and Unrau p. The whole discussion was plagued by the fact that often 
the relevant definitions and conventions are not stated explicitely. Unluckily the real 
or apparent differences left often the impression that the higher twist corrections are 
kind of undefined. This is not true. OPE gives a unique result and by now all groups 
working on this problem have reached the same results. In this paper we present our 
calculation in detail with all definitions given, such that no misunderstandings should 
remain. 
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As indicated the aim of this paper is to identify higher twist operators, which contribute 
at the next-to leading twist level. We do not however attempt to calculate 0{as) 
corrections to the Wilson coefficients. 

The starting point of the calculation is the forward virtual Compton scattering ampli- 
tude: 

T^,{q,P,S)=i J d''^e"i<{PS\T{j^{^)MO))\PS) . (1) 

The target nucleon has four momentum (P^ = ) and spin S'^ (S*^ = — , S ■ 
P = 0). The virtual photon has four momentum q^^ , [q^ = —Q^ < 0). The nucleon state 
is normalized covariantly {PS\P'S') = 2P°(27r^)5^(P — P')6ss'- Using the analyticity 
properties of T^^^ 

^ [T^uiqo + ie) - T^u{qo - ie)] = -27tW^u (2) 
and equation (Q) we find the following form for the hadronic tensor: 

= ^/ d'C e^'-HPSl m),juiO)]_ \PS) . (3) 

We define commutators and anticommutators as {A, B}^ = AB + BA and [A, = 
AB — BA. The antisymmetric part of the hadronic tensor can be decomposed as follows: 

<. = i£;...A^ (^(5i(x,Q2)+ff2(x,Q^))5'^-52(x,QV^) , (4) 
with = p ■ q, the Bjorken variable x = Q'^ jlv and ^^"^^ = 1. 

We are only interested in g\ and g^^ so we consider only the antisymmetric part of T^^y. 

Performing a Cauchy integration one obtains the connection between and 
which scematicaly reads |^]: 

W^^x^ dx = ^Jdu; T^Au)u^~^-^ , (5) 

with uj = 1/x. Equation (^) holds for n = only if T^(a;) vanishes fast enough for 
a; ^ oo. General considerations based on Regge theory suggest that the part of 
that determines gi goes rapidly enough to zero, but the part, which determines 172 may 
not Q. In this publication we assume the validity of equation (|^) for all n € Nq. 

Starting from equation (||) one obtains a representation of as a series in u) in the 
kinematical domain where | uj |< 1, i.e.: 

m 

Where from equation (|5|) it follows that: 

£ dx = i&(r/'^ • (7) 
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The only Feynman amplitudes which contribute to up to l/Q^ order are shown 
in Figure 1 |5|. To proceed further we follow the method proposed by Shuryak and 
Vainsthein ||5|. In their formalism contribution of Figure 1 to equation (|l|) has the 
form: 

T;t(g,P,5) = -^(^,)5: J d'i {PSm^fqj (^^^-^7, + ^,-^^^) ^/|0)|P5). 

^ (8) 
The symbol ^(^;^) denotes antisymmetrization of the indices and i'. Through this 
paper we use the Schwinger notation in which |^) denotes a formal eigenvektor of 
the coordinate operator X^, i.e.: 

X,\C)=C^\C), (9) 
while is the momentum operator which satisfies 

[P^,P.L = ^gGf^u. (10) 

In the coordinate basis acts as the covariante derivative 

(^W) = ''^^(^ - + 9A,{cm - eo • (11) 

In this representation the propagator for massless quarks can be written as 

SiC,0) = iC\^\0), 5(0,0 = -(o|^|0. (12) 

In addition to derive equation (^) we have used the following identities: 

e-i^-«P^ = (P, + q,)e-"^-^ , e-"^<\0) = |0) . (13) 

To avoid misunderstanding in the following we will denote by round brackets, e.g. 
(0|...|0, the matrix elements in the formal coordinate space, to be distinguished 
from square brackets denoting matrix elements between physical hadronic states, e.g. 
{PS\ . . . \PS). The Dirac spinor of the quark of flavor / is denoted ipf and qf is its 
charge. The sum runs over the light flavors (u,d,s) in the nucleon under consideration. 
Note that P denotes the impuls operator in these formulas, while q and ^ are c-numbers. 
The matrix element between proton states in (|8|) is taken at the scale Q^. To evaluate 
equation (§) we expand (^) and — ^): 

1 1 1 ^1 1 ^1 ^1 1 ^1 «.l ^1 



r+4 



1 1 1 ^1 1 ^1 ^1 1 ^1 ^1 ^1 , 

^ -1 - 1^1 - l^l^l - l^l^l^l + ■■■ ■ (14) 



r-4 
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The l.h.s. is some non-local operator, while the r.h.s. gives local differential operators 
by identifying = = + gA^^ with the covariant derivative (see Eq. (|Tl|)). 
Inserting only the first term of the expansion into equation yields in the lowest 
order. 

= dH{ps\mfq){i^.4iu-iu4i,)^fmps) 



X 



= 2ie^,A<x%E(^^iV^/(0)'?/7VV'/(0)|PS) 

= ie,..A^^(2 5:g2^f), (15) 

f 

where we have used the identity 

TmTaTi' - 7i^7A7M = 2ie^j.A(T7^7'^ , (16) 
(7^ = i7'^7-^7^7^) and have introduced the matrix element of the local operator 

(PS|V5/(0)7'^7V/(0)|P5) = 2S^af . (17) 

For in the lowest order this gives: 



1 I 
{gi{x) + g2ix)) dx = -a(°) 







\2{x)dx = 0, (18) 



with a^^^ ~ ^/^/ • This is the well known result: 

'5i(x)dx = ^aW . (19) 

The corrections to this sum rule are obtained by calculating the next term contributing 
to denoted by (5r^. 

'^^A = -4 E / d%psmi,fq) {n4mii. - lAmh^) ^m\ps) . (20) 

Using the equations of motion we getQ: 

(21) 



^ due to the EOM / d*^(c|V/ ■ ■ • fV-/ |o) = = / d*C(C| V'/f • ■ • '^Z |o) 
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In the above formula we suppressed / d'^S,{C\''Pf • • • V'/|0)- Now we decompose the Lorentz 
tensor of rank three P"p/^'y°"^5 ^^^^ parts with spin three to zero: 

^al3a ^al3a ^ a/3 p2 a 5 _ J_ p2 /? 5 _ J_ p2 a 5 

^5 spinS^^'^S spin 2 -^glJ ^11 ^^9 1 1 7 7 

+le"/5-A^(5^^^<5^5 ^ (22) 

and equivalent for P°^P^^'^ . G\s is the dual QCD field tensor 

GxS = ■^^XSafBG'^'^ . (23) 
The operators O are defined as follows: 

O5 %n 3 = I {P'^P^r + P'^P"!^ + P^P^^r 

+ P" P°7^ + pl^ ^"^7" + P" P^7"^ 7^ 

- — + p-igP^^a + p2^°'^^/3A ^5 ^ (24) 

18 ^ ' 
Ot%.n2 = \(2P"pPr -P'^P'^l^ -P^P-T)!"" 

(2p2g"/3^^ _ p2^/5a^a _ p2^aa^/3^ ^5 _ (25) 

The operators O^p^^g and O^^p^ns are the same as O5 ^spin 3 and O'^'^'^^^^^^ they just 
go without the 7^ matrix. The spin one part of equation (22) is checked to be the 
correct one by taking suitable traces, the spin zero part by contraction with £ai3(TX'- 
The spin three and two operators are traceless. The spin three part is known to have 
the form of equation Q. O5 ^s'pin 2 determined unique, because it is the only 

operator that fulfills the identity (|22|). Although this spin two operator has no specific 
symmetry properties, it can be written as a linear combination of two tensors with 
mixed symmetry. This corresponds to the fact, that a tensor of rank three contains two 
irreducible representations with spin two (as three vectors can be coupled to spin two 
in two different ways). 

Inserting equations (^) and ( p2[ ) into equation (^) gives: 



5T^ 



-{PS\^fmQj [ o"pl% + <n% + ^^g'-'P'r - -^^a'^'-P'!' 



E 

/ 



18" ' 6 

{PS\^f{o)q} (of;,, 3 + of:,,n 2 + ^a-'P'ri' - ^a'-'^PW 



(26) 
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For symmetry reasons only the following terms contribute to 5T^^ 



(27) 



Note that due to the contraction only O5 , but not O contribute to gi and g2 ■ Additional 
only that part of Of 2 ^^at is symmetric in a/? (Of ^j^^ 2) contributes. With 



O 



aficr 

5 spin 2 



2 ((1 + PQ/3)0f spin 2 + (1 - Pa/3) of spin 2 



afia 



5 spin 2 2 



+ ^((l-Pa/3)05 spin2y = 



(28) 



we get 



'^5 spin 2 



1 |^2P" P^7'^ + 2P'^ P"7'^ - P'^ P°''jl^ - P" ^^^7" - P'^7" - P" P'^7^^ 7^ 



-- (2p25°/37- - p2/-7" - p25"-7'3) 7^ . 



(29) 



is the operator, which changes the indices a and fi. As stressed by Shuryak and 
Vainshtein every step of this calculation is gauge invariant (0 = g^T^jy = q'^T^j^y). 

The matrix elements of the relevant local operators are as follows: 

{ps\4,f{Q)ofi^,i^fmps) = 



2a 



(2) 
/ 



M2 



- (P'^P'^S'' + P'^P'^S" + p'^p^s'^ + pPp'^S'' + P'^P'^S'^ + P-^P^^' 
6 ^ 

(g"^5'" + /'^5" + g"'^5^)] , (30) 



(P5|Vi/(0)Of;in2V'/(0)|P5) = 

= ^(P5|V;/(0) [7"5G'^'' + 7^5(5"''] V/(0)|PS) - traces 



= 2df 



1-0/(0) k^^G^'' + 7^5(5"''] V/(0)|PS) - traces = 

^ f Q^pOL pP gcr _|_ ^^pP pa ga _ pf3 pa ga pa pa gfl pa pa q(3 

.6 ^ 



/\/f2 , 



p<rpPga^ 

(31) 
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{PS\^f{0)P^-/^ro^PfmPS) = {PS\i>f{G)gG^p^^ijfmPS) = 2M^ffs^ . (32) 
With this definitions ^ for the matrix elements and equation (^) one gets from and 

J\x {g,{x,Q')+g,{x,Q')) = ia^ + ^ (a(2) + 4d(2) + 4/(^)) , 
/ dx g2{x,Q^) = 0, 

JO 

dx (51 (x, Q2) + 52(x, Q')) x' = ^a(2) + 1^(2) , 

f'dxg,{x,Q')x' = -ia(2) + >), (33) 
Jo 6 6 



with 

a(") = qjaf^ , = ^ qjdf^ , /(") = J] qjff"^ . (34) 

/ / / 

Including the next term of the expansion in would result in corrections of the order 
(Q^/M^)^ for the first moments, corrections of the order Q^/M^ for the third moments, 
and it would give values for the fifth moments. We can therefore conclude: 

l'dxg,{x,Q^) = \-^'^ + ^,{a^'^+^d('^+^f^'^)+Oi^), (35) 

dxg2{x,Q^) = + O(-^), (36) 

1 1 /if2 

dx5i(x,gV = 2''^'^+0{q2)^ (37) 
j\xg,{x,Q')x' = -\<^^'^ + \d^'^+0{^). (38) 



Equation ( |3q ) states that up to twist four and in the leading order of as there are 
no terms which contribute to the first moment of g2, i.e. which would violate the 
Burkhardt-Cottingham sum rule The sign in front of /^^^ in equation ( |35|) is 

positive. This is in agreement with the result of Shuryak and Vainshtein Q and Balitsky, 
Braun, Kolesnichenko [Q]. In their revised paper Ji and Unrau have derived the /^^^ 
term with a negative sign, but this is only due to their convention for D^, see footnote 
to equation (|3^). We want to stress that all three calculations agree with each other 
when the differences due to different conventions are taken into account. 

With a|f^^„t„n = 0.022 ± 0.002 from the E130 and the EMC [|] measurement, a'^l 

and a"^^ from the EMC measurement, and the estimates dproton ~ — 1.4- 10~^±4.0- 10~^ 

^Note the sign difFerence in equation (^) with respect to Ji and Unrau j^, owing to the fact, that 
they use = 9p + iff^p. 
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(if 



-0.026 ± 0.004) and /^^^ 



-0.050 ± 0.018 (/P 



- -u.u.iu _L u.uu-*; aiiKx Jproton " AT " 0-018 ± 0.018), wllich we 

get from the QCD sum rule calculation of Balitsky, Braun and Kolesnichenko 0] one 
obtaines 







(39) 

GeV2 



as(<5^)x . 33-8nfaJQ2) 

0.124 ■ (1 - ' ) + 0.013 • (1 ^ M - (0.018 ± 0.0072 

vr 33 — 2nf vr 



The QCD radiative corrections are taken from [|l3 
of the first moment of gi on for low (nj 
result of ref. Note, however, that in their publications Jbbk 



Figure 2 shows this dependence 
: 3). Our result agrees with the 
-7^ and therefore 



Sa BBK = ^Io-Ibbk-^ ~ The different sign in front of the corrections of 

equation ( ^9|) compared to the sign of Ji and Unrau (equation (29)) is due to fact, 

that although their calculation agrees with ours, numerically the bag model prediction 

(2) 

for the matrix element /proton opposite sign than the predictions from QCD sum 
rules. 
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Figure Caption 



Figure 1: Graphs contributing to T{jn{$,) 

Figure 2: Dependence of the first moment of on (solid line). The dotted line 
shows only the radiative corrections. The dashed line corresponds to — 0.018/Q^. 
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